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Thomas has recently derived scaUng laws for X-ray radiation from electrons accelerated in plasma 
bubbles, as well as a threshold for the self-injection of background electrons into the bubble [1 . To 
obtain this threshold, the equations of motion for a test electron are studied within the frame of the 
bubble model, where the bubble is described by prescribed electromagnetic fields and has a perfectly 
spherical shape. The author affirms that any elliptical trajectory of the form x'"^ /^^ -\- y''^ — is 
solution of the equations of motion (in the bubble frame), within the approximation p'y^ /p'x^ <^ 1. 
In addition, he highlights that his result is different from the work of Kostyukov et al. [2 , and 
explains the error committed by Kostyukov-Nerush-Pukhov-Seredov (KNPS). 

In this comment, we show that numerically integrated trajectories, based on the same equations 
than the analytical work of Thomas, lead to a completely different result for the self-injection 
threshold, the result published by KNPS [2 . We explain why the analytical analysis of Thomas 
fails and we provide a discussion based on numerical simulations which show exactly where the 
difference arises. We also show that the arguments of Thomas concerning the error of KNPS do 
not hold, and that their analysis is mathematically correct. Finally, we emphasize that if the KNPS 
threshold is found not to be verified in PIC (Particle In Cell) simulations or experiments, it is due 
to a deficiency of the model itself, and not to an error in the mathematical derivation. 



Authors of Ref. [T and Ref. [2^ have considered a 
model in which the bubble is described by prescribed 
electromagnetic fields and has a perfectly spherical shape 
in the laboratory frame, whose radius is and velocity is 

Vp = c^l — l/^p. They obtained different thresholds for 

electron self-injection into the bubble. Whereas Thomas 
argues that an error has been committed in the work of 
KNPS, leading to wrong conclusions, we will show in this 
comment that the conclusions of KNPS are correct (in 
the frame of the considered model) and that the math- 
ematical derivation of Thomas is erroneous. We begin 
by demonstrating that there is no elliptical solution for 
the equations of motion, whatever the initial conditions. 
Then, we explain why the arguments of Thomas concern- 
ing the error of KNPS do not hold, and we present numer- 
ical results showing agreement with the KNPS threshold. 
Finally, we provide a discussion based on numerical sim- 
ulations which show exactly why considering the trajec- 
tory as elliptical leads to erroneous conclusions. We give 
qualitative arguments which highlight that the consid- 
ered model could be too simple to quantitatively describe 
the self-injection physics. 



In the following, we use the prime to indicate quantities 
defined in the bubble rest frame, as opposed to quantities 
defined in the laboratory frame. In addition, quantities 
are normalized by the choice mg = c = e = ujp = l where 
ujp is the plasma frequency. Derivatives with respect to 
the electron proper time r are indicated with a dot: A = 
dA/dr. 



ELLIPTICAL TRAJECTORY 

In our conventions, the system of equations given by 
Eqs. (12) and (15) of Ref. [1 (equations of motion in 
the bubble frame) is written 



(1) 

(2) 



From these equations, Eq. (18) of Ref. [T can be estab- 
lished (with a minus sign instead of a plus sign in the 
l.h.s, and a factor rrie in the r.h.s) and is written 



1 d 



X -y 



IpV' dr 
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Note that, while Thomas made use of the approximation 
P'y^ /p'x^ <C 1 to derive Eq. O), this last equation can be 



derived without this approximation, such that, according 
to him, elliptical trajectories are not only approximate 
solutions (in the sense p'y^ /p'^^ <C 1) but exact solutions 
to the equations of motion. 

If Eqs. ([l]) and ([2| imply Eq. ([3|, the reverse is false. 
Providing initial conditions (x'(0), ^'(0), x'(0), ^'(0)) are 
known, there are an infinite number of solutions for Eq. 
(|3|, while only one for the system ([T]) + ([2|. Thomas states 
that ^^This equation is satisfied by any trajectory of the 
form x'^ ^ y''^ = R^^\ There is an infinite number 
of elliptical trajectories of this type, and they can be 
parametrized by x'ij) = jpR cos 0^ u' {^) = RsvciO where 
6> is a function of r (specifying a particular solution) . Ac- 
cording to the Thomas' affirmation, any elliptical trajec- 
tory is solution of Eq. (|3|, which means in mathematical 
terms: V6> G (x'(r) = -fpR cos 0,y\r) = RsinO) G 



2 



S'(3): where S(^s) is the solution space of Eq. ([3|. Insert- 
ing this parametrization into Eq. (|3| shows that terms 
in 6'^ and in do not cancel out, and that any trajectory 
of the form + y''^ = R'^ is not solution of Eq. (l3|. 

Instead, we obtain a differential equation for ^, which has 
an unique solution Os{r) providing that the initial con- 
ditions (6>(0),^(0)) are known. We note x^^ = jpRcosOs 
and y^ = Rs'mOs. By derivation, the trajectory (x^^, y^) is 
solution of Eq. ^ and satisfies x'^ j^^ -^y'^ = R? . How- 
ever, because Eq. ([3| is not equivalent to Eqs. (l) + (2), 
{^'s -til's) is Gi priori not a solution of the system (l) + (2). 
To show that {x'^^y'^) is effectively not a solution of the 
system Q+Q, its expression can be inserted in Eqs. ([T]) 
and ([2| . Here we propose a simpler demonstration based 
on a Taylor expansion of the solution of Eqs. ([T]) + ([2| 
around the initial time r = 0: 



x' = x\0) + x\0)t + -Y^^ + 6^ 



(4) 



(5) 



An elliptical trajectory has to satisfy the following rela- 
tion: 



xx' + jpyy 



0, 



(6) 



at all time r. The initial conditions, compatible with 
Eq. (§, are x'{0) = 0, y'(0) = n, i'(0) = p^o, 2/'(0) = 
0, where p'^Q is the only free parameter. To derive the 
second, third and fourth derivative of x' and y' at t = 
from Eqs. (jlj and the following relation is useful (it 
can be obtained from Eqs. (8) and (12) of Ref. [Tj): 



d 



-(7' + x')x'. 



We obtain from Eqs. Q and ([2]), and using Eq. 0: 
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(8) 



where 7o^ = We can insert this Taylor expansion 

of the solution of the equations of motion ( pi) a nd ([2| in 
the relation for the elliptical trajectory, Eq. (|6|), to check 
if the real solution is elliptical or not, in the limit r <C 1. 



FIG. 1. Electron trajectory for 7^ = 10 and for the initial 
conditions x'(0) = 0, y' {^) = 4, x'(0) = p'^Q, y{0) = [where 
Pxo is chosen as the solution of Eq. The numerical 

solution of Eqs. ([T]) and ^ is in solid red line, while the 
elliptical trajectory is in dashed green line. 



Identifying each order of expansion gives: 



Order : 

Order 1 : 
Order 2 : 

Order 3 : 



= 0, 



PxO 

= 0, 



v3^2 



0, 



(9) 
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ho+p'xo)p'xO 



il'o+P'xo) 



0. 



(10) 



Eqs. ([9| and (10) have to be satisfied simultaneously, 
whereas there is only one free parameter p^Q. These equa- 
tions are in fact incompatibles, they can not be satisfied 
simultaneously. For example, if we consider the limit 

-IpV^ from Eq. (g. 



(7) \Pxo\ > 1, we obtain p'^ 



xO 



which is not solution of Eq. (10) 



In addition to this analytical analysis, a numerical in- 
tegration of the equations of motion can be performed 
to verify if the trajectory can be elliptical, providing the 
correct choice of initial conditions. The value of p'^Q is 
chosen as the solution of Eq. (|9|, so that the trajectory 
is effectively elliptical to the lowest order in r. A nu- 
merically integrated trajectory is displayed on Fig. [l]for 
the parameters r5 = 4 and 7^ = 10. It is easily seen that 
the real trajectory does not follow an ellipse. We checked 
that errors due to finite time step and numerical trunca- 
tion were negligible; varying the time step or the level of 
truncation has no effect on the result. The trajectory is 
also found not to be sensitive to initial conditions, for the 
time scale of interest. Moreover, we performed a cross- 
verification of both the analytical and numerical calcula- 
tions. The Taylor expansion is valid only for r <C 1, and 
the time required for the electron to reach the back of the 
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FIG. 2. Electron trajectory for 7^ = 10 and for the initial 
conditions x{0) = 0, y{0) = 0.01, x (0) = p^o, ^ (0) = 
[where p^o is chosen as the solution of Eq. ^] . The numerical 
solution of Eqs. ([T]) and ^ is in solid red line, while the 
Taylor expansion of the solution, given by Eqs. Q, (jsj and 
([8|, is in dashed green line (superposed to the red linejT The 
elliptical trajectory is in dotted blue line. Trajectories are 
plotted up to r = 0.5. 

1 /3 

bubble is, in orders of magnitude, r ~ JpTb/lPxol^ ^ • 

1 /3 

Therefore we need <C 1 for the expansion to be valid 
on the length scale of interest (the bubble extension). On 
Fig. [2] is represented both the analytical Taylor expan- 
sion [given by Eqs. Q, ([5| and ([8|] and the numerically 
integrated trajectory [solution of Eqs. ([T]) and ([2|], up to 
r = 0.5, for n = 0.01 and 7^ = 10. The choice n = 0.01 
is only used here to perform a verification between an- 
alytical and numerical calculations (but this case does 
not have any physical relevance since the bubble model 
makes sense only for ao > 2, i.e. for > 2^/2). Both 
trajectories are very close to each other, confirming both 
the analytical and numerical calculations. 

We conclude that the real trajectory is not elliptical, 
whatever the initial conditions. We will see in the Discus- 
sion section why incorrectly considering the trajectory as 
elliptical leads to erroneous conclusions. 



ON THE ERROR OF KNPS 

Thomas argues that in the work of KNPS, the approx- 
imations made in Eqs. (4) to (7) of Ref. [2 are too 
restrictive and fail to correctly predict the self-injection 
threshold. This can be easily understood by regarding 
at Eq. (6) of Ref. [2 : X necessarily decreases, even 
when Px ^ 00 (in the notation of Ref. [2], X = ^/r5 = 
{x — Vpt)/ri)^ Px = Px/rl). Such equations can not de- 
scribe the injection, since when the electron is injected, 
X is increasing. Nevertheless, Eqs. (4) to (7) of Ref. 
are only used to obtain the numerical coefficient Px — 1.1 



at the moment where Py = for the first time, and to 
insert it into Eq. (3). The KNPS threshold is based on 
the conservation of the hamiltonian H between the ini- 
tial time and the critical time where Py = for the first 
time, and no approximation is needed in this approach. 
From that, Eq. (3) of Ref. [2] is established. Moreover, a 
simple analysis in orders of magnitude of the equations of 
motion Eqs. (1) and (2) of Ref. [2 shows that Px oc r^. 
Inserting this behavior in Eq. (3) of Ref. [2] demon- 
strates the self-injection threshold of Eq. (9) in Ref. [2], 
but without the numerical coefficient. The numerical co- 
efficient can then be evaluated by drastically simplifying 
the equations of motion, as done by KNPS. In reality, 
the coefficient could have a very weak dependance on the 
parameters r^j and 7^, since the real equations depend on 
them. 

Therefore, the arguments of Thomas concerning the er- 
ror of KNPS do not hold, and the semi-analytical deriva- 
tion of KNPS is correct. 



NUMERICAL THRESHOLD 

In order to verify the mathematical validity of the 
self-injection thresholds proposed either by Thomas or 
KNPS, we have integrated the equations of motion and 
scanned all the parameter space (r5,7p), with the same 
initial conditions as Thomas or KNPS, i.e. x'{0) = 0, 
y'{0) = r^, x'{0) = —jpVp^ y'{^) = (electron at rest in 
the laboratory frame). The electron is considered to be 
injected if = ^'^/7p + ^ ^5 at all time steps. Note 
that if the electron escapes the bubble before 7/ = 0, it 
will never come back inside, so that imposing the condi- 
tion r ^ r5 only when ^ ^ gives the same result. The 
numerical result is presented on Fig. [3] and is in agree- 
ment with the work of KNPS. In the frame of the model 
considered by KNPS and Thomas (with initial conditions 
corresponding to an electron at rest in the laboratory 
frame), the threshold is written > 1.307^. 

DISCUSSION 

In this comment, we have analyzed the mathematical 
validity of the results proposed either by Thomas and 
KNPS, considering a particular model where the bubble 
is considered perfectly spherical and described by pre- 
scribed electromagnetic fields. However, it is clear that 
such a simple model can potentially fail to correctly de- 
scribe the physical mechanisms present in the blow-out 
regime of laser-plasma interaction. In the work of KNPS 
[2 , only one PIC simulation has been performed, while 
several simulations with very different parameters should 
be performed to confirm the linear self-injection thresh- 
old. In addition, we highlight that, while Thomas consid- 
ers the parameter 7^ as a direct function of the electron 
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FIG. 3. For each value of rt and 7^, a value of 1 is displayed 
if the electron is injected and otherwise. The initial con- 
ditions are x\0) = 0, y\0) = x\0) = -jpVp, y\0) = 
(electron at rest in the laboratory frame), and a trajectory is 
considered injected if r ^ rb at all time steps. The frontier 
between injected and non-injected trajectories follows a line 
of equation jp — 0.77 rt (plotted in red line), whose numer- 
ical coefficient is very close to the value obtained by KNPS 
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FIG. 4. Numerically integrated electron trajectory for 7^ = 
200 and for the initial conditions x' (0) — 0, ^^(0) = 12, 
X {0) = —'jpVp, y{0) = (electron at rest in the laboratory 
frame), in solid red line. The ellipse of equation r = is in 
dashed green line, and the ellipse of equation r — ^max IS m 
dotted blue line. The inset displays a zoom of the back of the 
bubble and shows the point where = for the first time. 
At that point, r :^ rmax ^ rt and |x^|/7p ^ r^. 



plasma density 7^ ex , it should be instead con- 

sidered as the bubble back gamma factor which can be 
much lower due to the time evolution of bubble. Indeed, 
Kostyukov and co-workers included the rate of bubble ex- 
pansion in the bubble back gamma factor [3] and found 
similar results that those of Kalmykov et al. [T, which 
have explicitly studied electron injection inside a time- 
dependent bubble. This bubble back gamma factor has 
to be properly taken into account if we want to verify 
the linear self-injection threshold of KNPS by PIC sim- 
ulations. 

We have seen that, contrary to the Thomas' affirma- 
tion, there is no elliptical solution for Eqs. ([T]) and ([2|. 
The difference between the result of Thomas and the lin- 
ear self-injection threshold of KNPS can be understood 
as follows. In the bubble frame, the hamiltonian is writ- 
ten %' = ^' — (j)' ^ where (j)' 2jp(j) = — 7pr^/4, (j)^ and (j) 
being the scalar potential respectively in the bubble and 
laboratory frame. The model is time-independant in the 
bubble frame, therefore is conserved and 7'+7pr^/4 = 
7o + 7^^5/4 [this is Eq. (9) of Ref. [1 J. Because 
7' ^ 1, there is a maximal value for r, which is given 
by rmax ^ ^5(1 + 4/rg)^/2 for 7p > 1 and 7^ ^ 7^. For 
large values of r5, rmax becomes very close to r5, and 
if a small error is committed, an electron can be seen 
injected while it is not in the frame of the considered 
model. Thomas used Eq. ([7|), considering 7' as a con- 
stant (which is equivalent as saying that the trajectory 
is elliptical, according to the conservation of H')^ and 
studied the motion in terms of the and p'^ variables. 
In reality, 7' is not constant along the trajectory, which 



induces some degree of error in the calculation of the rela- 
tion between and p^, given by Eq. (20) of Ref. [1 . In 
fact, the difference in the result of Thomas arises when 
he considered an electron to be injected if x' ^ —Jp^b 
when = 0, implicitly assuming an elliptical trajectory, 
for which y' = {) and \x'\/^p = r when = 0, such 
that the condition {x^l/jp ^ r^ is equivalent to r ^ r^. 
But because the trajectory is not elliptical, when p'^ = 0, 
y' ^ and \x'\/jp 7^ r such that even if r^ < r ^ rmax 
(the electron is not injected), we can have Ix'l/jp ^ r^. 
For large r^, even if r = rmax when = 0, because rmax 
is very close to r^, \x'\/jp will be smaller than due to 
the non-zero value of y\ For example, for = 12 and 
7p = 200, according to Thomas the electron is injected, 
while it is not according to KNPS. Figure [4] displays the 
corresponding trajectory (with initial conditions for an 
electron at rest in the laboratory frame) and the ellipses 
of equation r = and r = rmax- At the moment where 
= for the first time, r rmax > r^ (it is considered 
as non-injected by KNPS), whereas Ix'l/jp ^ r^ (it is in- 
jected according to the Thomas' criterion). This example 
highlights that because r^ and rmax are very close to each 
other, a small error in the derivation or in the criterion 
can considerably change the conclusion (injected or non- 
injected). In addition, in that case, it is clear that 7' is 
not constant at all, since it almost attains 7' = 1 (when 
r rmax) and it attains very large values 7' 7q during 
the period where the electron is inside the bubble. 

We have performed a complete scan of the parameter 
space (r5,7p), as for Fig. [sj but applying either the con- 
dition r ^ r5 or \x'\/jp ^ r^ at the moment where p'^ = 
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FIG. 5. For each value of Vh and 7p, a value of 1 is displayed 
if \x'\/^p ^ Th at the moment where = for the first 
time, and otherwise. The initial conditions are x' (0) — 0, 
^^(0) = rb, x^(0) = —'jpVp, y{0) = (electron at rest in the 
laboratory frame). 



the KNPS threshold. But there is no physical meaning 
for applying \x'\/jp ^ as a criterion for self-injection. 
As we can see on Fig. |4j the electron streams backwards 
after the point where p'^ = and Ix'l/jp ^ r5. It should 
not be considered as injected in the frame of the con- 
sidered model. Therefore, the threshold for self-injection 
indicated in Eq. (22) of Ref. [1 is incorrect, because it 
relies on the elliptical trajectory which is in contradiction 
with the equations used to derive the threshold. 

Nevertheless, the present discussion emphasizes that, 
because when increases r^ax becomes very close to 
r^, a small deformation of the bubble structure, or the 
consideration of the field enhancement at the back of 
the bubble due to electron crossing, or the considera- 
tion of self-consistent screened fields, could considerably 
change the conclusion about injection or non-injection in 
the bubble. Considering these effects in the model and 
confirming or invalidating the KNPS result are areas for 
future works. 



for the first time. For the condition r ^ r^, the result is 
similar to Fig. |3] but with a slightly different numeri- 
cal coefficient, the threshold being > 0.957^. For the 
Thomas' condition, \x'\/jp ^ r5, the result is non-trivial 
and is displayed on Fig. [5] According to this criterion, 
self-injection occurs for much larger values of jp than for 
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